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Abstract. Using the light- front quark model, we calculate the transition form 
factors, decay rates, and longitudinal lepton polarization asymmetries for the exclusive 
rare B s — > (K, r]^)(£ + £~, vti>t {£ — e,fi,r) decays within the standard model, 
taking into account the r\ — rf mixing angle. For the mixing angle 9 — —20° 
(—10°) in the octet-singlet basis, we obtain BR(P S — > riJ2 v i^e) = 1-1 (1-7) x 10~ 6 , 
BR(B S -> 77/1+ = 1.5 (2.4) x 10~ 7 , BR(B S -> Tyr+r") = 3.8 (5.8) x 10~ 8 , 
BR(B S -> rjY, v lvi) = !-8 (1-3) x 10~ 6 , BR(B S -> rfn + n~) = 2.4 (1.8) x 10~ 7 , 
and BR(P S — > t)'t + t~) = 3.4 (2.6) x 10 -8 , respectively. The branching ratios for 
the B s — > K(v£V tl t + t~) decays are at least an order of magnitude smaller than those 
for the B s -> r/W {vm , £+l~) decays. The averaged values of the lepton polarization 
asymmetries for B s -> (K,rj^)£ + £- are obtained as (Pf) M = (P^ = (p/) M = 
-0.98, (Pf > T = -0.24, (P£) T = -0.20 and (p/) r = -0.14, respectively. 
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1. Introduction 

The study of the exclusive decays in the beauty sector allows one to explore the standard 
model (SM) and search for new physics effects. The B factory experiments such as BaBar 
at SLAC, Belle at KEK, LHCb at CERN, and B-TeV at Fermilab make precision tests 
of the SM and beyond the SM ever more promising. Especially, the 5 s -meson system 
becomes a key element in the 5-physics program of B factories ever since the first 
evidence for B s production at the T(55) was found by the CLEO collaboration [HE]. 
The DO [[3] and CDF [I] Collaborations have made measurements of the charge-parity 
(CP) violating weak B s — B s mixing phase <fi s in B s —y J/ip<p decays. While the SM 
expectation [5],[6] is nearly zero, the measured <p s differs from by more than 3cr (but 
with a sizable error). This measurement of cj) s inconsistent with zero (if confirmed) would 
indicate an evidence of new physics. Recently, the Belle Collaboration also measured 
the branching ratios of the B s —y J/ip<p an d B s — y J/i/jt) decays and the preliminary 
result [7] of the B s — y J/ip<p decay is about 3 times larger than that for the B s —y J/iprj 
decay. This ratio agrees with a rough estimate obtained within the naive quark model 
(neglect octet-singlet mixing), where the ss part of the rj meson wave function is one 
third in contrast to the fully ss content of mesons. With the upcoming chances that a 
numerous number of B s mesons will be produced at hadron colliders, one might explore 
the exclusive rare B s decays to (K, rj, r]')£ + £~ (and V(Vi) {£ = e,//, r) induced by the 
flavor- changing neutral current (FCNC) transitions b —y (d, s). Since in the SM the rare 
B s decays are forbidden at tree level and occur at the lowest order only through one-loop 
penguin diagrams [U [9j [101 [TTJ [121 HSl [111 [15] , the rare B s decays are well suited to test 
the SM and detect new physics effects. While the experimental tests of exclusive decays 
are much easier than those of inclusive ones, the theoretical understanding of exlcusive 
decays is complicated mainly due to the nonperturbative hadronic form factors entered 
in the long distance nonperturbative contributions. Therefore, a reliable estimate of the 
hadronic form factors for the exclusive rare B s decays is very important to make correct 
predictions within and beyond the SM. The r\ — rf mixing angle may also be extracted 
from the rare B s decays to rj and rf final states. 

In our previous work [IB] , we have analyzed the exclusive rare B —y K£ + £~ 
decays within the framework of the SM, using our light-front quark model (LFQM) 
based on the QCD- motivated effective LF Hamiltonian [T71 [T8l [T9] . The experimental 
values of the branching ratios BR(5 -> K£ + £r) = (0.75^;! ± 0.09) x 10" 6 from 
Belle [20] and (0.34 ± 0.07 ± 0.02) x 10~ 6 (£ = e,/x) from BABAR [21] detectors are 
consistent with our LFQM prediction 0.5 x 10~ 6 [16] based on the SM. Recently, we 
also analyzed B c properties and various exclusive decay modes such as the semileptonic 
B c —y (D,r] c , B, B s )£ve decays [22], the rare B c — y D^t r £~ [23], and the nonleptonic 
two-body B c — y (-D( s ), r) c , B^)(P, V) decays |24j (here P and V denote pseudoscalar 
and vector mesons, respectively). The form factors f±(q 2 ) and friQ 2 ) for the exclusive 
rare decays [23] between two pseudoscalar mesons are obtained in the Drell-Yan-West 
(q + = q° + q 3 = 0) frame [25] (i.e., q 2 = — < 0), which is useful because only 
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the valence contributions are needed unless the zero-mode contribution exists. The 
covariance (i.e., frame independence) of our model has been checked by performing the 
LF calculation in the q + = frame in parallel with the manifestly covariant calculation 
using the exactly solvable covariant fermion field theory model in (3 + 1) dimensions. We 
also found the zero-mode contribution to the form factor f~{q 2 ) and identified [22J the 
zero-mode operator that is convoluted with the initial and final state LF wave functions. 

The purpose of this paper is to extend our our LFQM [T6|fT7tfT8|fT9l[22|[23l[2i]to 
calculate the hadronic form factors, decay rates and the longitudinal lepton polarization 
asymmetries (LPAs) for the exclusive rare B s — > (K, rj, i]')£ + £~ and V(pi decays 
within the SM. The LPA, as another parity-violating observable, is an important 
asymmetry [26] and could be measured at hadron colliders such as LHCb. In particular, 
the r channel would be more accessible experimentally than e- or //-channels since the 
LPAs in the SM are known to be proportional to the lepton mass. There are some 
theoretical approaches to the calculations of the exclusive rare B s — > r/£ + £~ [27J [23 [29] 
and B s — > rj'£ + £~ [2HI [29] decays, but not the B s — > K£ + £~ decay mode as far as we 
know. 

The paper is organized as follows. In Sec. 2, the SM operator basis, describing 
the b — > (d, s)(£ + £~, vgpi) transitions, is presented. In Sec. 3, we briefly describe the 
formulation of our LFQM and the procedure of fixing the model parameters using the 
variational principle for the QCD motivated effective Hamiltonian. We discuss the rare 
decays between two pseudoscalar mesons using an exactly solvable model based on the 
covariant Bethe-Salpeter (BS) model of (3 + l)-dimensional fermion field theory and 
show the equivalence between the results obtained by the manifestly covariant method 
and the LF method in the q + = frame. We then present the LF covariant forms of the 
form factors f±(q 2 ) and friQ 2 ) obtained from our LFQM. The n — rj' mixing angle for the 
B s — > r]V> transitions is also discussed in this section. In Sec. 4, our numerical results, 

1. e. the form factors, decay rates, and the LPAs for the rare B s — > (K,r),r)')(£ + £~ , Vfpi) 
decays are presented. Summary and discussion of our main results follow in Sec. 5. 

2. Effective Hamiltonian 

In the SM, the exclusive rare B s — > P q (£ + £~, v^pi) (q = d,s) decays are at the quark 
level described by the loop b — > q (£ + £~ , i>£Ve) transitions, and receive contributions from 
the Z(7)-penguin and 1^-box diagrams as shown in Fig. [TJ 

The effective Hamiltonian responsible for the b — > q£ + £~~ (q = d, s) decay processes 
can be represented in terms of the Wilson coefficients, Cj S , C| ff and C w as (9] 



where Gf is the Fermi constant, a em is the fine structure constant, and Vy are the 
Cabibbo-Kobayashi-Maskawa (CKM) matrix elements. The relevant Wilson coefficients 
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Figure 1. Loop diagrams for B s — > P q (£ + £ , vivg)(q — d, s) transitions. 

Ci can be found in Ref. [9]. The effective Hamiltonian responsible for the b — > qv^vi (q = 
d, s) decay processes is given by [TU [12] 

Kf = V^W* - 75)^(1 - TB>t, (2) 

l\J 2n sin u\y 

where x t = {m t /M w ) 2 and X(x t ) is the top quark loop function [HI [12], which is given 
by 

, w . x /2 + x 3x — 6 , \ 

X ( x )=ol 7 + 7 T^ lnx • ( 3 ) 



,2 — 1 (x — 1 

Besides the short distance (SD) contributions, the main effect on the decay comes from 
the long distance (LD) contributions due to the cc resonance states (J/i[),i/;', ■ ■ •)• The 
effective Wilson coefficient C| ff taking into account both the SD and LD contributions 
has the following form [9] 

Cf(s) = C 9 + Y SD {s) + Y LD {s), (4) 

where the explicit forms of Y S d(s) and Y LD (s) can be found in [HI ED]- For the LD 
contribution Y LD (s), we include two cc resonant states J/ip{lS) and ip'{2S) and use 
r(J/> -»■ = 5.26 x 1CT 6 GeV, M Jfi , = 3.1 GeV, r w = 87 x 1CT 6 GeV for 

J/ip(lS) and rr>' -»■ ^) = 2.12 x 1CT 6 GeV, = 3.69 GeV, I> = 277 x 1CT 6 
GeV for ip'(2S) [31]. 

The LD contributions to the exclusive B s — > P q (q = d, s) decays are contained 
in the meson matrix elements of the bilinear quark currents appearing in and 
T-L v c q 1 ■ In the matrix elements of the hadronic currents for B s — > P q transitions, the 
parts containing 75 do not contribute. Considering Lorentz and parity invariances, 
these matrix elements can be parametrized in terms of hadronic form factors as follows: 

J» ee (P q \rfb\B.) = f + (q 2 )P» + /_(?V, (5) 

and 

J£ = (P q \qia^q„b\B s ) = J^f^ [fP* ~ - Mj> )g"], (6) 

where P = Pb s + Pp q and q = Pb s — Pp q is the four-momentum transfer to the lepton 
pair and 4mj < q 2 < (M Bs — M Pq ) 2 . We use the convention = (i/2)[y fl , j u ] for the 



Exclusive rare B s — y (K, rj, n')£ + £ decays in the light-front quark model 5 

antisymmetric tensor. Sometimes it is useful to express Eq. (jSJ) in terms of f+(q 2 ) and 
fo(q 2 ), which are related to the exchange of 1~ and + , respectively, and satisfy the 
following relations: 

,2 



/+(o) = / (o), fo(g 2 ) = U(q 2 ) + 



q 



Ml - Ml 



f-(q 2 )- 



(7) 



With the help of the effective Hamiltonian in Eq. (PQ) and Eqs. (E]) and (JHJ), the transition 
amplitude M = (P q £ + £-\H eS \B s ) for the B s -> 



M 
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decay can be written as 
The differential decay rate for B s — y P q £ + £~ is given by J32J [33] 
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m 2 /Mj> , m b = m h /M Bs and r = Mf, /M| . Equation <Q may 



|Cio| 
H = (s - 1 
with s = q 2 /M 2 Ba} t 

be written in terms of (/+, fo, fr) instead of (/+, /_, fr) as discussed in [16]. Note also 
from Eqs. (jUJ) and f flUj) that the form factor f~(q 2 ) does not contribute in the massless 
lepton limit. 

The differential decay rate for B s —y P q v$yt can be easily obtained from the 
corresponding formula Eq. fl9]) for B s —y P q i 
C10 — y X(x t )/ sin 2 6 W , i.e. 



by the replacement rhi — y 0, Of 1 —y 0, 



and Cf 



d 



MlG 2 



-0t o 



V tb V tq \ 2 <p¥\X(x t )\ 2 \f, 



11) 



ds^~ uiui 3-2 8 vr 5 sin 4 ^ 

where the factor of 3 in the numerator corresponds to the sum over the three neutrino 
flavors. Dividing Eqs. ([9]) and ( TTTT) by the total width of the B s meson, one can obtain 
the differential branching ratio dBR(B s -> P q £ + £~)/ds = (dT(B s -» P/ + £-)/r tot )/<is. 
As another interesting observable, the LPA, is defined as 

dTh=-i/ds — dTh=i/ds 



where h 



L ^ ^ dTh=-i/ds + dTh=i/ds' 
T(— 1) denotes right (left) handed £~ in the final state. From Eq. 



obtains for B s -y P q £ + £ 



Pl(s) 



2 1 — 4-! 



1/2 



UReCf - 



2m b CS 



H 1 + 2! + 



(12) 



one 



(13) 



Because of the experimental difficulties of studying the polarizations of each lepton 
depending on s and the Wilson coefficients, it would be better to eliminate the 
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dependence of the LPA on s, by considering the averaged form over the entire 
kinematical region. The averaged LPA is defined by 

3. Review of our LFQM 

The key idea in our LFQM [17, 18, 22J for the ground state mesons is to treat the radial 
wave function as a trial function for the variational principle to the QCD-motivated 
effective Hamiltonian saturating the Fock state expansion by the constituent quark and 
antiquark. The QCD-motivated effective Hamiltonian for a description of the ground 
state meson mass spectra is given by 



H qq = H + V qq = sjm? q + k 2 + sjm\ + k 2 + V qq , (15) 

where 

V qq = V + V hyp = a + br n -^ + ^^i V V coul . (16) 

or 3 m q m q 

In this work, we use the Coulomb plus linear confining (i.e. n = 1) potential together 
with the hyperfine interaction (S q ■ S q ) = 1/4 (—3/4) for the vector (pseudoscalar) 
meson, which enables us to analyze the meson mass spectra and various wave-function- 
related observables, such as decay constants, electromagnetic form factors of mesons in 
a spacelike region, and the weak form factors for the exclusive semileptonic and rare 
decays of pseudoscalar mesons in the timelike region [161 EH HI (El 1221 [231 1211 EH [35] . 

The momentum-space LF wave function of the ground state pseudoscalar mesons 
is given by \l/(xj, kjj_, Aj) = 7Z\ l \ 2 (x i , k i± )cf)(xi, kjj_), where 4>(x i ,'ki±) is the radial wave 
function and 7^AiA 2 is the covariant spin-orbit wave function. The model wave function 
is represented by the Lorentz-invariant variables, Xi = pf/P + , kjj_ = pjj_ — XjPj_ and 
Xi, where P" = (P + ,P-,P ± ) = (P° + P 3 , (M 2 + P 2 L )/P+,P ± ) is the momentum of 
the meson M, and pf and Aj are the momenta and the helicities of constituent quarks, 
respectively. 

The covariant form of the spin-orbit wave function for pseudoscalar mesons is given 

by 

y/2 ^M 2 -( mi -m 2 ) 2 

where Mq = S?=i(k 2 j_ + m 2 )/xi is the boost invariant meson mass square obtained from 
the free energies of the constituents in mesons. For the radial wave function 
the Gaussian wave function: 

4-7T 3 / 4 fdk 
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Table 1. The constituent quark masses [GeV] and the Gaussian parameters /3[GeV] 
obtained from the linear potential in [18], which are necessary for B s — > (K, 77,77') 
decay modes, q — u and d. 



m q 


m s 


m b 




Psa 




0.22 


0.45 


5.2 


0.3886 


0.4128 


0.5712 



where (3 is the variational parameter. When the longitudinal component k z is defined 
by k z = (x — 1/2)Mq + [jn\ — iti\)/2Mq, the Jacobian of the variable transformation 
{x, k ± } -)• k = (k_L, k z ) is given by 



dk z M Q f [m 2 — m 2 " 12 



dx Ax 1X2 



Mo 



(19) 



The normalization factor in Eq. ({TBI) is obtained from the following normalization of the 
total wave function: 

1 dx [ ^p-\V(x,k i± )\ 2 = 1. (20) 



J 167T 3 

We apply our variational principle to the Q CD-motivated effective Hamiltonian first 
to evaluate the expectation value of the central Hamiltonian H + V , i.e., (<p\(H -\-Vo)\<p) , 
with a trial function (^(x^kjjj that depends on the variational parameter fl. Once the 
model parameters are fixed by minimizing the expectation value {4>\(Hq + Vo)|0), the 
mass eigenvalue of each meson is obtained as M qq - = (<p\(H + V qq -)\(p). Minimizing 
energies with respect to /3 and searching for a fit to the observed ground state meson 
spectra, our central potential Vq obtained from our optimized potential parameters 
(a = -0.72 GeV, b = 0.18 GeV 2 , and a s = 0.31) [17J for the Coulomb plus linear 
potential was found to be quite comparable with the quark potential model suggested 
by Scora and Isgur [36], where they obtained a = —0.81 GeV, b = 0.18 GeV 2 , and 
a s = 0.3 ~ 0.6 for the Coulomb plus linear confining potential. A more detailed 
procedure for determining the model parameters of light- and heavy-quark sectors can 
be found in our previous works ITS] . Our model parameters (m q ,(3 q q) obtained 
from the linear potential model relevant to this work are summarized in Table [TJ 
The predictions of the ground state meson mass spectra can be found in [22]. We 
should note that our model parameters (m, 0) automatically satisfies the normalization 
of the total wave function and were fixed by the variational principle to the QCD- 
motivated effective Hamiltonian. Those parameters in turn automatically satisfies the 
normalization of the electromagnetic form factors at q 2 = and every other physical 
observables obtained from our LFQM such as decay constants and electroweak form 
factors are the predictions. This distinguishes our LFQM from other quark model. 

3.1. Form factors for the rare B s — > P decays in our LFQM 



Most popular phenomenological LFQM uses the Gaussian wave function as the radial 
wave function due to its predictive power of various physical observables. However, 
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Figure 2. The covariant diagram (a) corresponds to the sum of the LF valence 
diagram (b) defined in < k+ < P+ region and the nonvalence diagram (c) defined in 
P 2 + < fc+ < region. The large white and black blobs at the meson-quark vertices 
in (b) and (c) represent the ordinary LF wave functions and the non-wave-function 
vertex, respectively. The small black box at the quark-gauge boson vertex indicates 
the insertion of the relevant Wilson operator. 



since the LFQM using the Gaussian wave function does not have a counterpart of 
manifestly covariant model, it is hard to check the covariance of the model. To check 
the covariance of LFQM, one can start from the manifestly covariant field theory 
model. For example, using the exactly solvable covariant BS model of (3+l)-dimensional 
fermion field theory [371 EH [39] , one can perform the LF calculation in parallel with 
the manifestly covariant calculation and compare the results from the two models. 
Comparing the LF results and the manifestly covariant results, we were able to derive the 
LF covariant form factors between two pseudoscalar meson explicitly and to analyze the 
zero-mode complication. Since the detailed procedure of finding LF covariant transition 
form factors (/+,/_, fa) was already given in our previous works [221 122], we shall 
briefly describe the essential procedure of obtaining the LF covariant form factors from 
the exactly solvable covariant BS model of (3+l)-dimensional fermion field theory and 
show the results of the LF covariant form factors. 

The covariant diagram shown in Fig. [21(a) is in general equivalent to the sum of the 
LF valence diagram 2(b) and the nonvalence diagram 2(c). The matrix element 
obtained from the covariant diagram of Fig. [2](a) is given by 

where g\ and #2 are the normalization factors which can be fixed by requiring both 
charge form factors of pseudoscalar mesons to be unity at zero momentum transfer, 
respectively. To regularize the covariant fermion triangle-loop in (3+1) dimensions, we 
replace the point gauge-boson vertex 7 M by a non-local smearing gauge-boson vertex 
(K x 2 /N Kl )^(K 2 2 /N k2 ), where N Al =p\- Ai 2 + ie and N A2 = p% - A 2 2 + ie, and thus 
the factor (A1A2) 2 appears in the normalization factor. Ai and A2 play the role of 
momentum cut-offs similar to the Pauli-Villars regularization [37J. Our replacement 
of 7 M by the non-local smearing gauge-boson vertex remedies the conceptual difficulty 
associated with the asymmetry appearing if the fermion loop were regulated by smearing 
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the qq bound-state vertex. The rest of the denominators in Eq. (I2ip . i.e., NiNj.N 2 , are 
coming from the intermediate fermion propagators in the triangle loop diagram and are 
given by 

N k = k 2 - m 2 + ie, Nj = p 2 - m 2 + ie {j = 1,2), (22) 

where mi, m, and m 2 are the masses of the constituents carrying the intermediate four- 
momenta V\ — P\ — k, k, and p 2 = P 2 — k, respectively. Furthermore, the trace terms 
S 11 from the vector current and St^ from the tensor current are given by 

S" = Tr[ 75 (^i + m 1 )7"(^ 2 + m 2 ) 75 (- f6 + m)] 

= A{pi (p 2 ■ k + m 2 m) + p 2 (p 1 ■ k + mini) + F(mim 2 - p\ ■ p 2 )}, (23) 

and 

S£ = Tr[7 5 (^ + m^io^ q v {i> 2 + m 2 )j 5 (- ft + m)] 

= - 4{p1[m(p 2 ■ q) + m 2 (k ■ q)] - pi 2 l [m(p 1 ■ q) + m 1 (k ■ q)] 

+ k**[m 1 (p 2 - q) -m 2 (pi- q)]}, (24) 

respectively. By doing the integration over k~ in Eq. (1211) . one finds the two LF time- 
ordered contributions to the residue calculations corresponding to the two poles in 
k~ , the LF valence contribution [Fig. EJ^b)] defined in < k + < P 2 region and the 
nonvalence contribution [Fig. [2(c)] defined in P 2 < k + < P^ region. The nonvalence 
contribution [Fig. Ws)] in the q + > frame corresponds to the zero mode (if it exists) 
in the q + — > limit [30]. Performing the LF calculation of Eq. f l2T|) in the q + = frame 
in parallel with the manifestly covariant calculation, we use the plus component of the 
currents to obtain the form factors f+{q 2 ) and friq 2 )- For the form factor f-{q 2 ), we 
use both the plus and perpendicular components of the currents. As we have shown 
in [22l [23], while the form factors f+{q 2 ) and fr{q 2 ) can be obtained only from the valence 
contribution in the q + = frame without encountering the zero-mode contribution, the 
form factor f-{q 2 ) receives the zero mode. In our recent analysis of semileptonic B c 
decays [22], we identified the zero-mode operator that is convoluted with the initial and 
final state LF valence wave functions to generate the zero-mode contribution to the form 
factor f-{q 2 ) in the q + = frame. Our method can also be realized effectively by the 
method presented by Jaus [39] using the orientation of the LF plane characterized by the 
invariant equation u-x = 0, where u is an arbitrary light-like four vector. More detailed 
analysis of the zero-mode operator and the LF covariance of the form factors f± and 
fx can be found in [221 E3]. While the manifestly covariant BS model of fermion field 
theory model is good for the qualitative analysis of the exclusive rare decays, it is still 
semi-realistic. We thus replace the LF vertex functions in the BS model with the more 
phenomenological Gaussian radial wave functions in our LFQM since the zero-mode 
operator is independent from the choice of radial wave function as discussed in [22] . 

The LF covariant form factors f±(q 2 ) and friq 2 ) for B s (qiq) — > P(q 2 q) transitions 
obtained from the q + = frame are given by (see [221 [22] for more detailed derivations) 

f + (q ) = dx {AiA 2 + k ± • kjj, (25) 

Jo J 167r3 Ij2 + k 2 ± l A 2 + k /2 
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f-(q 2 



1 f d 2 k ± 0i(x,k 

(1 — x)dx 



q ■ P 



12 



x(l - x)M{ 



k ± — m\vn + (m2 — m)^4i + 2 



k , +2 (k^ 



+ 2 



'k± • q±) 2 k_L • q± 



— [M 2 2 - (1 - x)(g 2 + g • P) + 2xM 2 



(1 - 2x)M 1 2 - 2 (mi - m)(mi + m 2 )] 



(26) 



/ T (g 2 ) = (Mx + M 2 ) ^(1 - x)dx J 



d 2 k± 0!(x,kj 



(x,k'_, 



167T 3 



A\ + k 2 



A 2 , + k' 2 



A\ - (mi - m 2 ) 



(27) 



where k'j_ = kj_ + (1 — x)q±, Ai — (1 — x)m,i + xm (i = 1,2), and q- P = M\ — M| with 
Mi and M 2 being the physical masses of the initial and final state mesons, respectively. 
Our results for the form factors given by Eqs. (!25|) - (!27|) are essentially the same as 
those presented in [41J. We should note that the LF covariant form factor f-(q 2 ) in 
Eq. ( 126]) is the sum of the valence contribution f™ l (q 2 ) and the zero-mode contribution 
f^ M '(q 2 ) [22]. Since the form factors f±(q 2 ) and /t(<? 2 ) are defined in the spacelike 
(q 2 = — < 0) region, we then analytically continue them to the timelike q 2 > region 
by changing q^ to — q 2 in the form factors. We also compare our analytic solutions with 
the double pole parametric form given by 

/<(0) 



1 - 0i s + a 2 s 2 
where 01 and o~ 2 are the fitted parameters. 



(28) 



3.2. r] — rj mixing for the B s — > rf'> decays 

In this subsection, we discuss the 77 — rj mixing to obtain the B s — > ry>'> transition form 
factors. The octet-singlet mixing angle 9 of rj and 77' is known to be in the range of —10° 
to —23° [31]. The physical 77 and 77' are the mixtures of the flavor SU(3) octet 77s and 
singlet ?7o states: 

/ 77 \ = / cos# -sin9 \ / r]s \ . v 

\rf J V sin ^ cos ^ / V Vo J ' 

where 77s = (uu + dd — 2ss)/\/Q and 770 = (uu + dd + ss)/\/3. Analogously, in terms of 
the quark-flavor(QF) basis 77^ = (uu + dd)/v2 and 77^ = ss, one obtains [12] 

/ 77 \ = ( cos0 -sin0 W r) q \ 
\ rj J \ sin cos / \ 77^ J ' 

The two schemes are equivalent to each other by (ft = 9 + arctany^ when SU/(3) 
symmetry is perfect. Although it was frequently assumed that the decay constants 
follow the same pattern of state mixing, the mixing properties of the decay constants 
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will generally be different from those of the meson state since the decay constants only 
probe the short- distance properties of the valence Fock states while the state mixing 
refers to the mixing of the overall wave function [42J. 

Defining (P(p)\J^\0) = —ifpp^ (P = r], 7/) in the QF basis, the four parameters 
fp and fp can be expressed in terms of two mixing angles (<f> q and S ) and two decay 
constants (f q and / s ), i.e. [42J, 

n n \ = ( 008 ^ - sin ^ \( /« M m 

/} ) \ sin0 g cos0 s )\0 f s )■ 1 } 

The difference between the mixing angles g — S is due to the Okubo-Zweig-Iizuka(OZI)- 
violating effects [33] and is found to be small (<f) q — 4> s < 5°). The OZI rule implies that 
the difference between (fi q and S vanishes (i.e., 4> q = 4> s = 0) to leading order in the 
1/N C expansion. Similarly, the four parameters fp and fp in the octet-singlet basis 
may be written in terms of two angles (0 8 and 9 ) and two decay constants (/ 8 and / ). 
However, in this case, ^8 and turn out to differ considerably and become equal only 
in the SU/(3) symmetry limit [4~2| |4^]. 

We shall use the QF basis with the single mixing angle to analyze the B s — y rf'' 
decay modes. In this case, a generic form factor F and the branching ratio for the 
B s — y 7]^ are given by 

F B S ^ V { V ') = _ gin ^ CQS ^pBs^Vs ; (32) 

with the physical r/^ mass and 

BR[B S -> r](7]')£ + £~} = sin 2 0(cos 2 <j))BR[B s -> r] s £ + £-}, (33) 

respectively. Recently, the KLOE Collaboration [15] extracted the pseudoscalar mixing 
angle <fi in the QF basis by measuring the ratio BR(0 — y 7/7) /BR(0 — y 777). The 
measured values are = (39.7 ± 0.7)° and (41.5 ± 0.3 stat ± 0.7 syst ± 0.6 t h)° with and 
without the gluonium content for 77', respectively. However, since the mixing angle for 
i] — r\' is still a controversial issue, we use unspecified value for rather than adopting 
some specific value. 

4. Numerical results 

In our numerical calculations for the exclusive rare B s —y (K,7],r]')(i , £i7 £ ,£ + £~) decays, 
we use the model parameters (m q , 0) for the linear confining potential given in Table [U 
Although our predictions [22] of ground state heavy meson masses are overall in good 
agreement with the experimental values, we use the experimental meson masses [31] in 
the computations of the decay widths to reduce possible theoretical uncertainties. 

Note that in the numerical calculations we take (m c , m&) = (1.8,5.2) GeV in all 
formulas except in the Wilson coefficient Cg S , where (m c ,m^ po i e ) = (1.4,4.8) GeV have 
been commonly used [9] . For the numerical values of the Wilson coefficients, we use the 
results given by Ref. [9]: 

d = - 0.248, C 2 = 1.107, C 3 = 0.011, 
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q 2 [GeV 2 ] q 2 [GeV 2 ] 

Figure 3. The weak form factors for B s — > K(\eit panel) and B s —> r/ s (right panel) 
decays, respectively. 



Table 2. Results for form factors at q 2 = of B s — > (K,r],r]') transition and 
parameters Oi defined in Eq. (I28|) . The coefficients in rj and rf represent quark mixing 
angles, i.e. c v = — sine/) and c v ' = cos0, respectively. 



Mode 


Mo) 


o-i 




/-(0) 


O-l 


(T 2 




Mo) 




0-2 


B 3 


0.230 


-1.650 


0.822 


-0.201 


-1.638 


Oi 


135 


-0.228 


-1.633 


0.835 


B S ^T] 


0.291c,, 


-1.574 


0.751 


-0.231c^ 


-1.582 


Oi 


^25 


-0.280c^ 


-1.561 


0.782 




0.291cy 


-1.575 


0.770 


-0.225c^/ 


-1.570 


Oi 


135 


-0.300c v 


-1.561 


0.802 



C 4 = - 0.026, C 5 = 0.007, C 6 = -0.031, 

Cf = - 0.313, C 9 = 4.344, C w = -4.669, (34) 

and other input parameters are |X^ 6 X^*| = 0.039, |V^ 6 ^| = 0.008, = 129, 

M w = 80.43 GeV, m t = 171.3 GeV, sin 2 9 W = 0.2233, and t b o = 1.470 ps. 

In Fig. [31 we show the q 2 dependences of the form factors f±(q 2 ) and /t(<? 2 ) for 
the B s —7- K (left panel) and B s — > r] s with physical masses of 77 and 77' (right panel), 
respectively. The form factors at q 2 = and the parameters Cj of the double pole form 
defined in Eq. ( 128]) are listed in Table [2j The form factor f+(q 2 ) for the B s — > r] s has the 
same q 2 dependence (apart from the mixing angle 0) for both rj (solid line) and 7/ (circle) 
since / + does not depend on the daughter meson mass as one can see from Eq. ( 125|) . 
On the other hand, the form factors /_ and Jt between 77 and 77' are slightly different 
(apart from the mixing angle 0) since they involve the daughter meson mass as one can 
see from Eqs. ( 126]) and (l27j) . The form factors at the zero recoil point (i.e., q 2 = g^ax) 
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- B s -> K v V 




Figure 4. Differential branching ratios for B s — >• K Y vgvi (left panel) and B s — > 
rjs v lVt (right panel) decays, respectively. 



correspond to the overlap integral of the initial and final state meson wave functions. 
The maximum recoil point (i.e., q 2 = 0) corresponds to a final state meson recoiling 
with the maximum three- momentum \Pf\ — (M§ s — M 2 )/2Mb s in the rest frame of the 
B s meson. As a sensitivity check of our LFQM, we find for the B s —y K transition 
that our form factors are changed only about 3% as the light <i-quark mass varies about 
20%. This indicates that the transition form factors for b —y d decay processes are 
quite stable on the variation of d quark mass. The form factors for the B s —y n s have 
also been computed by Geng and Liu [28] using the similar LFQM but only with the 
valence contributions in the purely longitudinal q + ^ frame. Although the form factors 
/+ and /t at the maximum recoil point obtained from [28J do not receive nonvalence 
contributions, they receive nonvalence contributions at other nonzero q 2 values. The 
nonvalence contributions to /_ are more serious for the entire q 2 range including the 
q 2 = point. For instance, /_ obtained from [28] (see Fig. 2 in [28]) shows a sharp 
increasing as q 2 near the zero recoil point in contrast to our result. This indicates 
the nonvalence contribution to /_ is quite large, which in particular overestimate the 
branching ratio for the r dilepton decay mode. Although the form factor f-(q 2 ) does 
not contribute to the branching ratio in the massless lepton (£ = e or fi) decay, it is 
important for the heavy r decay process. 

We show our results for the differential branching ratios for B s —y (K, r] s ) J2 v i^t 
with physical masses of 77 and 7/ in Fig. H] and B s —y (K,rj,r]')£ + £^ (£ = fx and r) 
in Fig. |5l respectively. For the B s —y (K, rj, rj')£ + £~ transitions in Fig. EJ the solid 
(dashed) lines represent the results without (with) the LD contribution to Cg S . Since 
the B s — y K is induced by b —y d transition compared to the B s — y r) s induced by 
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0.9 



0.9 



0.7 



Figure 5. Differential branching ratios for B s — > K£ + £~ (upper panel) and B s 
i]£ + £" (middle panel) and B s — > r]'£ + £~ (lower panel) with £ — fi and r, respectively. 
The solid and dashed lines represent the results without and with the LD contributions, 
respectively. 
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b — > s at the quark level, the branching ratios for the final K meson are order of 
magnitude smaller than the corresponding branching ratios for the final 77 meson. For 
the B s — > (K,r],r]')£ + £~ decays (see Fig. [5]), the LD contributions (dashed lines) clearly 
overwhelm the nonresonant branching ratios near J/%p(lS) and ip'(2S) peaks, however, 
suitable £ + £~ invariant mass cuts can separate the LD contribution from the SD one away 
from these peaks. This divides the spectrum into two distinct regions [26], H6]: (i) low- 
dilepton mass, Am 2 < q 2 < M 2 ,^ — 5, and (ii) high-dilepton mass, M 2 , + 5 < q 2 < g^ax; 
where 5 is to be matched to an experimental cut. 

Our predictions for the nonresonant branching ratios are summarized in Table [3] 
with general form of the mixing angle <ft in the QF basis. Our results are also compared 
with other theoretical predictions such as the LF and constituent QM [28] and the 
QCD sum rules (SR) [29J within the SM. Since the amplitude B s — )■ (K, 77, r]')£ + £~ is 
regular at q 2 = 0, the transitions B s — > (K,r],r]')e + e~ and B s — > (K,r],rj')fi + /i~ have 
almost the same decay rates, i.e. insensitive to the mass of the light lepton. Our 
predictions of branching ratios are close to the QCD SR results [29] but a bit smaller 
than the LFQM results [28]. But the results from [28J could be lowered if the nonvalence 
contributions are properly taken into account. For the mixing angle 9 = —20° (—10°) 
in the octet-singlet basis, which corresponds to <fi = 34.74° (44.74°) in the QF basis, we 
obtain BR(B S ryE^) = 1-1 (1-7) x KT 6 , BR(B S ->• ri/i+fi-) = 1.5 (2.4) x 1(T 7 , 
BR(B S r/T+T-) = 3.8 (5.8) x KT 8 , BR(B S 7/ £ v e j?t) = 1.8 (1.3) x 10" 6 , 
BR(B S ///1+/1-) = 2.4 (1.8) x 10- 7 , and BR(B S r/'r+r-) = 3.4 (2.6) x KT 8 , 
respectively. 

It is also worth comparing the branching ratios between B s — >■ K and B s — > 77, 
which may be written as 

BR(B S -» Kfi+fi-) _ 1 V td 2 

BR(B S ^ ¥ V1 " ^0 Ws {1 - Asu ^ (35) 
where the SU (3) correction term A S u(3) is estimated about 0.3 in our model calculation. 
Such a kind of relation may be further scrutinized by considering an additional correction 
term neglected in the effective Hamiltonian as discussed in [29J. The branching ratios 
with the LD contributions for B s — > (K, 77, r]')£ + £~ (£ = /n, r) are also presented in 
Table H] for low- and high-dilepton mass regions of q 2 . 

In Fig. [6] we show the LPAs for B — > (K,rj,r]')£ + £~ {£ = /x, r) as a function of s. 
In both fi and r dilepton decays, the LPAs become zero at the end point regions of s. 
However, we note that if me = 0, the LPAs are not zero at the end points. As in the 
case of the B ->- K^^T pH [33J, [32l [47] and B c ->• D {s) n + ^~ [23] decays where P L ~ -1 
away from the end point regions, the LPAs away from the end point regions are also 
close to —1 for the B s — > (K,i],i]')fi + transitions. In fact, the Pi for the muon decay 
is insensitive to the form factors, e.g. our Pi away from the end point regions is well 
approximated by 
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Table 3. Nonresonant branching ratios (in units of 10 7 ) for B s — > (K 1 ri,rj')'^2 l UiVi 
and B s — > (ri,ri')£ + g~ transitions compared with other theoretical model predictions 
within the SM. 



Mode 


This work 


m 


[29] 


D i IS V* j/„f/„ 

d s — t I\ 2^ VlVl 


1 m 








OO. 1 bill G3 


Oo.o bin <p^±jr^ivij 


^Zl .0 33 41.0 J bill (Jy^bei A J 






04.1 bill </)(L^MJ 


(ou.l ± lo. yj bin </)(bet d) 




9K 9 nna 2 ^1 
L\J.L COb (/J 




^lU.U 31 O.U J COb </J ^bet rv J 






oy. ( cob (/^ov^ivij 


^oo.y in o.y j cob </j ^beu d i 


O s 7 l\ jj, jjj 


fl 1 A 

U. l^i 






B s — > r)u, + u~ 


4.75 bin 2 


8.53 bin 2 0(LFQM) 


(2.73 ±0.68) bin 2 0( S et A) 






7.78 bin 2 0(CQM) 


(5.92 ± 1.59) bin 2 (bet B) 


B s ->■ t]'fi + ^ 


3.53 cob 2 


6.06 cob 2 0(LFQM) 


(1.96 ±0.53) cob 2 0(bet A) 






5.69 cob 2 0(CQM) 


(3.92 ± 1.07) cob 2 (bet B) 


B s -»■ Kr+r- 


0.03 






P s —7- r]T + T~ 


1.17bin 2 


1.67 bin 2 0(LFQM) 


(0.68 ±0.11) bin 2 0( S et A) 






1.67 bin 2 0(CQM) 


(1.82 ±0.34) bin 2 (bet B) 


B s -»■ r/V + r" 


0.51 cob 2 


0.83 cob 2 0(LFQM) 


(0.28 ±0.05) cob 2 0(bet A) 






0.72 cob 2 0(CQM) 


(0.69 ±0.13) cob 2 (bet B) 


Table 4. Branching i 


atios with the LD contributions for B c —> (D, D s )£ + £ for low 


and high dilepton mass regions of [GeV 2 ] 


obtained from the linear (HO) potential 


parameters. 






Mode 


Am} < q 2 < 8.5 


14.5 < q 2 < gLx 


B s -> Kfi + fi~ 


7.72 (6.63) x 10" 9 


2.27 (2.62) x 10" 9 


B s -)• Kr+r 






2.43 (2.66) x 10~ 9 




2.86 (2.44) 


bin 2 X 10 7 


7.00 (8.10) bin 2 x 10~ 8 


P s —7- 7]T + T^ 






9.13 (9.85) bin 2 x 10~ 8 


B s -)■ T]'fi + fl- 


2.54 (2.17) 


COb 2 x 10~ 7 


2.54 (2.13) cob 2 x 10~ 8 


B s -> r/r+r- 






3.42 (3.67) cob 2 x 10" 8 



in the limit of C| ff —> from Eq. (fl3|) . It albo bhowb that the Pl for the \i dilepton 
channel ib inbenbitive to the little variation of C% s as expected. On the other hand, the 
LPA for the r dilepton channel ib benbitive to the form factorb. Similar obbervation hab 
albo been made in our recent work for B c — > (D, D s )£ + £~ decayb [23J. 

The averaged valueb (Pf'' 7< ' ) )^ of the LPAb for B s -> (K,r]^)£ + £r without the LD 
contributions are (Pf )„ = (P2)„ = (P/) M = -0.98, (Pf ) T = -0.24, (P^) T = -0.20 
and (P2 ) t = —0.14, rebpectively. 
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Figure 6. Longitudinal lepton polarization asymmetries for B s — > K£ + £~ (upper 
panel), B s — > i]£ + £~ (middle panel) and B s —> r)'£ + £~ (lower panel). The same line 
codes are used as in Fig. [5j 
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5. Summary and Discussion 

In this work, we investigated the exclusive rare semileptonic B s — y (K,r),r) ; )(i , ei7e,£ + £~) 
(£ = e, fi, t) decays within the SM, using our LFQM constrained by the variational 
principle for the QCD motivated effective Hamiltonian with the linear plus Coulomb 
interaction [TTt ITS] . Our model parameters obtained from the variational principle 
uniquely determine the physical quantities related to the above processes. This approach 
can establish the broader applicability of our LFQM to the wider range of hadronic 
phenomena. The weak form factors f±(q 2 ) and fril 2 ) for the B s —y (K,r],r]') decays 
are obtained in the q + = frame (q 2 = — < 0) and then analytically continued to 
the timelike region by changing q^ to — q 2 in the form factors. The covariance (i.e., 
frame independence) of our model has been checked by performing the LF calculation 
in parallel with the manifestly covariant calculation using the exactly solvable covariant 
fermion field theory model in (3 + l)-dimensions. While the form factors f+{q 2 ) and 
/r(<? 2 ) are immune to the zero modes, the form factor f-{q 2 ) is not free from the zero 
mode. Using the solutions of the weak form factors obtained from the q + = frame, 
we calculated the branching ratios for B s —y (K,T),7)')(v£Ve,£ + £~) and the LPAs for 
B s — y (K, rj, i]')£ + £~ including both short- and long-distance contributions from the 
QCD Wilson coefficients. Our numerical results for the nonresonant branching ratios of 
B s —y rj^'^viVi, t + t~) decays are C(10~ 6 , 1CT 7 , 1CT 8 ) in orders of magnitude, 

respectively. The branching ratios for the B s —y K(v£Ue,£ + £~) decays are at least an 
order of magnitude smaller than those for the B s — y rf^iytf/^ £ + £~) decays. The averaged 
values (P^ )^ of the LPAs for B s — y (K, r}^)£~^t~ without the LD contributions are 
(P?)r = (PI), = (Pi), = "0.98, (Pf > T = -0.24, (P£) T = -0.20 and (P/) T = -0.14, 
respectively. These polarization asymmetries provide valuable information on the flavor 
changing loop effects in the SM. Of particular interest, we estimated that the ratio 
BR(P S — y K[i + fi~)/BR(B s — >■ r]fi + fi~) differs from the SU/(3) symmetry limit (apart 
from the mixing angle) by about 30%. Such a kind of relation may help in determining 
the rj — T)' mixing angle. 

Acknowledgments 

This work was supported by the Korea Research Foundation Grant funded by the Korean 
Government(KRF-2008-521-C00077). 

References 

[1] Artuso M et al (CLEO Collaboration) 2005 Phys. Rev. Lett. 95 261801 
[2] Bonvicini G et al (CLEO Collaboration) 2006 Phys. Rev. Lett. 96 022002 

[3] Abazov V M et al (DO Collaboration) 2007 Phys. Rev. Lett. 98 121801; 2008 Phys. Rev. Lett. 
101 241801 

[4] Aaltonen T et al (CDF Collaboration) 2008 Phys. Rev. Lett. 100 161802 
[5] Lenz A and Nierste 2007 J. High Energy Phys. 06 072 
[6] Bona M et al 2006 J. High Energy Phys. 10 081 



Exclusive rare B s — > (K, 7], t]')£ + £ decays in the light-front quark model 



19 



[7] Drutskoy A 2009 larXiv:0905. 29591 

[8] Grinstein B, Wise M B, and Savage M J 1989 Nucl. Phys. B 319 271 

[9] Buras A J and Miinz M 1995 Phys. Rev. D 52 186 

[10] Misiak M 1993 Nucl. Phys. B 393 23; ibid. 439, 461(E) (1995) 

[11] Inami T and Lim C S 1981 Prog. Theor. Phys. 65 297 

[12] Buchalla G, Buras A J, and Lautcnbacher M E 1996 Rev. Mod. Phys. 68 1125 

[13] Ali A, Manned T, and Morozumi T 1991 Phys. Lett. B 273 505; Ali A Acta Phys. Pol. B 27, 3529 
(1996) 

[14] Kim C S, Morozumi T, and Sanda A I 1997 Phys. Rev. D 56 7240 

[15] Aliev T M, Kim C S, and Savci M 1998 Phys. Lett. B 441 410 

[16] Choi H M, Ji C R, and Kisslinger L S 2002 Phys. Rev. D 65 074032 

[17] Choi H M and Ji C R 1999 Phys. Rev. D 59 074015 

[18] Choi H M and Ji C R 1999 Phys. Lett. B 460 461; 1999 Phys. Rev. D 59 034001 

[19] Ji C R and Choi H M 2001 Phys. Lett. B 513 330. 

[20] Abe K et al (Belle Collaboration) 2002 Phys. Rev. Lett. 88 021801 

[21] Aubert B et al (BABAR Collaboration) 2006 Phys. Rev. D 73 092001 

[22] Choi H M and Ji C R 2009 Phys. Rev. D 80 054016 

[23] Choi H M 2010 Phys. Rev. D 81 054003 

[24] Choi H M and Ji C R 2009 Phys. Rev. D 80 114003 

[25] Drcll S D and Yan T M 1970 Phys. Rev. Lett. 24 181; West G 1970 Phys. Rev. Lett. 24 1206 

[26] Hewett J 1996 Phys. Rev. D 53 4964; Kruger F and Sehgal L M 1996 Phys. Lett. B 380 199 

[27] Skands P Z 2001 J. High Energy Phys. 01 008 

[28] Ceng C Q and Liu C C 2003 J. Phys. G 29 1103 

[29] Carlucci M V, Colangelo P, and De Fazio F 2009 Phys. Rev. D 80 055023 

[30] Faessler A, Gutsche Th, Ivanov M A, Korncr J G, and Lyubovitskij V E 2002 Eur. Phys. J. direct 
(74, 18 

[31] Amsler C et al (Particle Data Group) 2008 Phys. Lett. B 667 1 

[32] Melikhov D, Nikitin N, and Simula S 1998 Phys. Rev. D 57 6814 

[33] Ceng C Q and Kao C P 1996 Phys. Rev. D 54 5636 

[34] Choi H M 2007 Phys. Rev. D 75 073016 

[35] Choi H M 2008 Phys. Rev. D 77 097301 

[36] Scora D and Isgur N 1995 Phys. Rev. D 52 2783 

[37] Bakker B L G, Choi H M, and Ji C R 2001 Phys. Rev. D 63 074014; 2002 Phys. Rev. D 65 

116001; 2003 Phys. Rev. D 67 113007. 

[38] de Melo J P B C and Frederico T 1997 Phys. Rev. C 55 2043; de Melo J P B C, Frederico T, 

Pace E, and Salme G 2006 Phys. Rev. D 73 074013 

[39] Jaus W 1999 Phys. Rev. D 60 054026 

[40] Choi H M and Ji C R 1998 Phys. Rev. D 58 071901 (R); 2005 Phys. Rev. D 72 013004; Brodsky 
S J and Hwang D S 1999 Nucl. Phys. B 543 239; M. Burkardt 1993 Phys. Rev. D 47 4628; de 
Melo J P B C, Sales J H O, Frederico T, and Sauer P U 1998 Nucl. Phys. A 631 574c 

[41] Cheng H Y, Chua C K and Hwang C W 2004 Phys. Rev. D 69 074025 

[42] Fcldmann T, Kroll P, and Stcch B 1998 Phys. Rev. D 58 114006; 1999 Phys. Lett. B 449 339 

[43] Feldmann T 2000 Lnt. J. Mod. Phys. A 15 159 

[44] Leutwyler H 1998 Nucl. Phys. B(Proc. Suppl.) 64 223 

[45] Ambrosino F et al (KLOE Collaboration) 2007 Phys. Lett. B 648 267 

[46] Ali A, Guidice G F, and Mannel T 1995 Z. Phys. C 67 417 

[47] Roberts W 1996 Phys. Rev. D 54 863; Burdman G 1995 Phys. Rev. D 52 6400 



